Magnetization of Graphane by Dehydrogenation 
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Using first principles calculations, we show that each hydrogen vacancy created at graphane 
surface results in a local unpaired spin. For domains of hydrogen vacancies the situation is, however 
complex and depends on the size and geometry of domains, as well as whether the domains are 
single- or double-sided. In single-sided domains, hydrogen atoms at the other side are relocated to 
pair the spins of adjacent carbon atoms by forming 7r-bonds. Owing to the different characters of 
exchange coupling in different ranges and interplay between unpaired spin and the binding geometry 
of hydrogen, vacancy domains can attain sizable net magnetic moments. 
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Graphene[l|, a truly two-dimensional (2D) crystal of 
honeycomb structure, has sparked considerable interest 
not only because of its charge carriers behaving like mass- 
less Dirac fermions[2|,|3|,|4j], but also the unusual magnetic 
roperties displayed by its flakes and nanoribbons0, 0, 
H, [^, [13] • In addition to numerous experimental and 
theoretical studies on the physical properties of graphene, 
efforts have been also devoted to synthesize various types 
of derivatives of graphene. More recently, a 2D hydro- 
carbon material in the family of honeycomb structure, 
namely graphane is synthesized [11]. Interesting prop- 
erties such as reversible hydrogenation-dehydrogenation 
with changing temperature ll| , the electronic structure 
with a wide band gap [13, | have been revealed soon af- 
ter its synthesis. In this letter, we reveal that graphane 
can be magnetized by dehydrogenation of domains on 
its surfaces. Large magnetic moments can be attained 
in a small domain on the graphane sheet, depending on 
whether the defect region is one-sided or two-sided. Our 
predictions are obtained from the state-of-the-art s pin 
polarized first- principles plane- wave calculations (Til. Il5j ] 
within the LDA noncollinear calculations including spin- 
orbit interaction, using (llxllxl) supercells. Details of 
our method can be found in [l6|, [l7| • 

Graphane, in its chair conformation as illustrated in 
Fig. [TJa), is derived by the adsorption of a single hydro- 
gen atom to each carbon atom alternating between the 
top (A) and bottom (B) side in the honeycomb struc- 
ture. A charge of 0.1 electrons is transferred from H to C 
leaving behind positively charged H atoms on both sides 
of a double layer of negatively charged (-0.1 electrons) 
C atoms. Graphane having a 2D quadruple structure 
has the work function $=4.97 eV, which is ~ 0.2 eV 
larger than that of graphene. In contrast to semimetal- 
lic graphene, graphane is a semiconductor with a wide 
direct band gap of 3.42 eV calculated by LDA but cor- 
rected to be 5.97 eV with GWo self-energy method, as 
shown in Fig. [TJb). Doubly degenerate states at the T- 
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point at the top of the valence band are derived from 
2p x - and 2pj,-orbitals of carbon atoms. The edge of the 
conduction band is composed mainly from G-p z orbitals. 
Calculated phonon bands all having positive frequencies 
confirm the stability of 2D graphane. High frequency 
vibration modes associated with C-H bonds are well sep- 
arated from the rest of the spectrum, in Fig. [He). 

The creation of a single H-vacancy at the hydrogen 
covered surfaces gives rise to the spin polarization in 
the non-magnetic perfect graphane. Desorption of a sin- 
gle H atom from graphane is an endothermic reaction 
with 4.79 eV energy. Various techniques, such as laser 
beam resonating with surface- hydrogen bond[l8|], strip- 
ping with ionic vapor [19| and scission of C-H bonds with 
subnanometer Pt clusters [201] . can be used to create H- 
vacancy(ies). Upon desorption of a single hydrogen atom, 
local bonding through sp 3 hybrid orbital is retransformed 
into planar sp 2 and perpendicular p z (tt) orbitals. At the 
vacancy site one unpaired electron accommodated by the 
dangling p z orbital contributes to the magnetization by 
one \ib (i-e. Bohr magneton). The exchange interaction 
between two H- vacancies calculated in a (llxllxl) super- 
cell is found to be non-magnetic for the first and second 
nearest neighbor distances due to spin pairings. Since the 
7T-7T interaction vanishes for farther distances, antiferro- 
magnetic (AFM) state between two H-vacancies for the 
third and fourth nearest neigbor distance is energetically 
favorable. The occurrence of long range spin interactions 
in carbon based structures was explained before by the 
superexchange 21 and magnetic tail interaction [22l]. 

As for the islands of H-vacancies at the single (top) side 
of graphane, we consider various geometrical domains, 
where H atoms at their edges and inside are removed as 
seen in FigJ^l For a triangular domain specified as A| at 
the top side, H atoms attached to three carbon atoms lo- 
cated at each edge are removed. Hydrogen atom which is 
normally adsorbed on the central C atom at the bottom 
side moves to the corner. Under these circumstances, 
spins of three hydrogen-free C atoms are antiferromag- 
netically ordered to yield a net magnetic moment of 1 
fiB- Noncollinear calculations with spin-orbit interaction 
fix the directions of spins, which are tilted relative to 
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FIG. 1: (Color online) (a) Top and side views of atomic struc- 
ture showing of graphane primitive unit cell with Bravais lat- 
tice vectors bi and b2 and buckling of alternating carbon 
atoms, A and B, in honeycomb structure S, bond lengths 
dc-c and dc-H optimized using LDA. Large green (light) 
and small orange (dark) balls indicate C and H atoms, re- 
spectively, (b) Energy band structure is calculated by using 
LDA and corrected using GWo (shown by blue lines and or- 
ange dots). For graphene, linear band crossing at Dirac point 
is shown by dashed grey lines, (c) Calculated phonon bands 
and density of states DOS projected to C and H atoms. 



the normal to the graphane plane. For A|, a triangu- 
lar domain has ten H atoms removed from the top side 
of graphane. While part of six H atoms are attached to 
carbon atoms from bottom are relocated, remaining two 
H atoms are released by forming H2 molecule. At the 
end spins are paired and the net magnetic moment of the 
domain becomes vanished. Generally, for a small single- 
sided domain, if N t , the total number of H atoms 
stripped, is an even number so that adjacent 7r-orbitals 
form spin paired 7r-bonds. In this case, H atoms below 
the domain are relocated (without facing any energy bar- 
rier) to pair adjacent 7r-orbitals to form maximum num- 
ber of 7r-bonds. At the end, a large buckled regions in- 
side the domain tends to be flattened and reconstructed 
to make nonmagnetic graphene-like planar structure. In 
the case of A§, while spins are paired through the for- 
mation of 7r-bonding between two adjacent C atoms fol- 
lowing the relocation H atoms at the bottom side, the 
unpaired spins at the corner atoms are aligned in the 
same direction to yield a net magnetic moment of //t=3 
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FIG. 2: (Color online) Calculated magnetic state of various 
domains of single-sided H-vacancies, where all H atoms at- 
tached to C atoms from upper side in the unshaded region (de- 
lineated by dash-dotted lines) including edges, are removed. 
The triangles are specified by A° with n indicating the max- 
imum number of C atoms at one edge and s signifies the 
single-sided dehydrogenation. Similar symbols are used also 
for hexagonal, ff| and lane L s n (n =4,5) domains. Total mag- 
netic moment \it and its components fi x , [i y and [i z are given 
in units of the Bohr magneton [ib- Magnetic moments on 
C atoms are shown by red (black) arrows. Relocations of 
H atoms at the other side of graphane are shown by curly 
arrows. For the sake of clarity 7r-bonds formed after the relo- 
cation of bottom H atoms are indicated only for A J, LJ and 
LI structures. 



fiB- The tendency to pair the spins of adjacent C atoms 
to form 7r-bonds are seen better in lane domains. Let 
us consider L\ and L§ in FigO Because of relocation of 
H atoms at the bottom side, two pairs of nearest neigh- 
bor C atoms form 7r-bonds and hence pair their spins. 
At the end, L% has /zt=0. For L§ having odd number 
of H-vacancy, while two pairs of C atoms are bound by 
two 7r-bonds, C atom at the center has an unpaired spin 
and attains /iT=l Ms- In a similar manner, the hexago- 
nal domain i/f has total of seven C atoms at its center 
and corners, all H atoms stripped from top side. At the 
bottom side, H atoms are relocated and hence the spins 
of adjacent C atoms are paired to result in a total net 
magnetic moment of n=l \ib- 

We next show in FigJ3] that the magnetic moment of 
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FIG. 3: (Color online) Net magnetic moments in Bohr magne- 
ton within the triangular A„, hexagonal H„, rectangular R„ 
and lane L„ domains, which are delineated by dash-dotted 
lines and have n carbon atoms at their edges. Here d signifies 
the double-sided dehydrogenation. Random shaped domain 
including both one and two-sided H-vacancy parts is also il- 
lustrated. 



graphane can be tuned by changing the size and geom- 
etry of a given double-sided H-vacancy domain. In this 
case the situation is not complex and allows us to figure 
out the magnetic moment of the entire structure easily. 
Based on noncollinear calculations including the spin- 
orbit coupling, the direction of the unpaired spins on the 
^4-type C atoms freed from H atoms is found to be op- 
posite to that of the spins of 5-type C atoms. However, 
instead of AFM spin ordering, lowest energy state of lane 
defects consisting of even number of C atoms is NM due 
to the entirely paired p z orbitals. Also, large double-sided 



domains including lane defects with equal number of A- 
and -B-type C atoms are found to be NM. The result- 
ing net magnetic moment of a double-sided H-vacancy 
domains can be given by /it = (N t — Ni,)hb, where N t 
and Nb denote the number of stripped H atoms from the 
top and bottom sides, respectively. Accordingly, the net 
magnetic moment induced in Af, A3, Af and Ag do- 
mains are 2, 3, 4 and 5 the [Ib respectively. The same 
argument can be applied to rectangular R%, hexagonal 
and lane domains. Even the magnetic moment of 
a domain having arbitrary shape including various single- 
sided and double-sided H-vacancy parts, as indicated in 
FiglH can be retrieved by the arguments discussed above. 
Non-integer value of /it is due to severe distortion of 
structure. We also note that our results regarding to the 
unpaired spin of a domain and their net magnetic mo- 
ment are in compliance with Lieb's theorem[23j, which 
distinguishes A- and i?-sublattices in honeycomb struc- 
ture. 

In conclusion, we showed that the interaction between 
unpaired spins associated with H vacancies in graphane 
gives rise to interesting magnetic structures. We revealed 
simple physical mechanisms underlying the magnetism 
of single-sided and double-sided vacancy domains. For 
single-sided domains, owing to the tendency to pair the 
spins of 7r-orbitals of adjacent C atoms, some of the ad- 
sorbed H atoms at the bottom side are relocated. At 
the end, the net magnetic moments can be attained in 
vacancy domains depending on their size and shape. For 
double-sided domains, interactions underlying the gen- 
eration of net magnetic moment are relatively straight- 
forward and are in good agreement with Lieb's theorem. 
Since the exchange coupling between different domains 
are hindered by domain walls, very dense data storage 
can be achieved through uniform coverage of identical 
domains. It is also noted that a graphane flake compris- 
ing a domain with large magnetic moment can be utilized 
as a non-toxic marker for imaging purposes. While mag- 
netic 2D systems attract a great deal of attention due to 
their tunable properties at nanoscale, our results suggest 
that the size and ordering of magnetic moments of hy- 
drogen vacancy domains with thin walls can be used for 
future data storage and spintronics applications. 

Computing resources used in this work were partly 
provided by the National Center for High Performance 
Computing of Turkey (UYBHM) under grant number 2- 
024-2007. 
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